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Auger recombination rates in nitrides from first principles 
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We report Auger recombination rates for wurtzite InGaN calculated from first principles density- 
functional and many-body-perturbation theory. Two different mechanisms are examined - inter- and 
intra-band recombination - that affect different parts of the emission spectrum. In the blue to green 
spectral region and at room temperature the Auger coefficient can be as large as 2xf0~'^''cm®s~^; in 
the infrared even larger. Since Auger recombination scales with the cubic power of the free-carrier 
concentration it becomes an important non-radiative loss mechanism at high current densities. Our 
results indicate that Auger recombination may be responsible for the loss of quantum efficiency that 
affects InGaN-based light emitters. 

PACS numbers: 71.20.Nr, 72.20.Jv, 79.20.Fv, 85.60.Bt 
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Indium gallium nitride (InGaN) alloys are now already 
being used for light emitting and laser diodes in the green 
to ultraviolet part of the spectrum^ but increases in in- 
ternal quantum efficiency (IQE) are still required to al- 
low broader applications. The IQE of InGaN devices is 
limited by loss mechanisms that, at high drive currents 
(i.e., high carrier concentrations) lead to a decrease in 
IQE, a phenomenon commonly referred to as "efficiency 
droop". The precise nature of these loss mechanisms has 
been the subject of intense debate, and a variety of can- 
didates have been proposed (see Ref.0). Recently, Shen 
et al. suggested Auger recombination as the dominant 
sourcei^i^ Loss due to Auger recombination scales with 
the cubic power of the free-carrier density and would thus 
dominate at the high carrier concentrations at which the 
reduction in IQE is observed. While Shen et al. found 
a cubic dependence of the IQE on the free-carrier con- 
centration in optically pumped InGaN LED devices,^ it 
is difficult to discriminate between different radiation- 
less processes experimentally. In this Letter we demon- 
strate by means of rigourous first-principles calculations, 
in which the Auger process can explicitly be isolated, that 
Auger recombination is indeed an important loss mecha- 
nism in wurtzite InGaN. 

In the direct Auger process, an electron recombines 
with a hole, but instead of emitting a photon the process 
results in the excitation of another carrier to a higher- 
energy state (see insets of Fig. [T|). This can also be 
viewed as two electrons colliding in the vicinity of a hole, 
resulting in a radiationless e-h recombination event, the 
energy and momentum of which is absorbed by the sec- 
ond electron (ee/i process). Alternatively, the hhe process 
involves two holes and one electron. A first-principles 
description of the Auger recombination rate therefore re- 
quires an accurate calculation of the band structure and 
the transition probabilities for all relevant eeh and hhe 
processes as inputs. 

We obtain the band structures of the host ma- 
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FIG. 1: eeh Auger coefficient for light holes as a function of 
fundamental band gap for a simulated InGaN alloy (see text) 
at T = 300 K and n = 1 x 10^^ cm"^. The statistical error bars 
of the Monte Carlo integration are smaller than the symbols 
for all data points presented. Intraband Auger, which involves 
scattering of an electron to an unoccupied part of the lowest 
conduction band (left inset, dashed line), dominates in In-rich 
alloys. However, for alloy concentrations relevant for solid- 
state lighting, interband Auger recombination, which involves 
scattering to the second-lowest conduction band (right inset, 
solid line), is dominant. 



terials by combining density-functional theory (DFT) 
with many-body perturbation theory in the GqWq 
approximation!^ This approach accurately describes 
band structures of solids as measured by direct 
and inverse photoemission/^i^ii^ Combined with the 
exact- exchange optimized effective potential approach 
(OEPx(cLDA)) for the DFT part, it produces accurate 
band structures for AIN, GaN and InN in their zinc- 
blende and wurtzite phasesi^ii^ We do not include spin- 
orbit coupling, which has only very small effects on the 
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band structure of nitrides. The results for wurtzite are 
shown in Fig. (5] 

The Auger rate R is given hy^^^^ 

S{ksuni)S{Esum) dkidk2dk3dk4 (1) 

where states 3 and 4 are, e.g. for eeh, electrons in the con- 
duction band and states 1 and 2 are a hole in the valence 
band and an electron in a higher-energy conduction-band 
state, ksum and E'sum are short for ki -t- k2 — ka — and 
El + E2 — E3 — E4, respectively. The statistics factor 

A,2,3,4 = (1 - / (El)) (1 - / {E2)) f (Es) f {Ei) (2) 

is determined by the Fermi occupation functions /, and 
-^1,2,3,4 are the Auger matrix elements of the screened 
Coulomb potential W 

Mi,2,3,4 = {c^ih\W\c^zcj)4) + EX, (3) 

where EX denotes exchange terms, as detailed in Ref . fisl. 

We compute Mi_2,3,4 directly with the OEPx(cLDA) 
wave functions and the non-local, dynamically screened 
Coulomb interaction Wq from an RPA calculation. 
Our OEPx(cLDA) calculations were performed with the 
S/PHI/nX plane-wave pseudopotential code,— while for 
the Go Wo calculations we have employed the gwst space- 
time codeJ^ii^iii With all operators being treated in real 
space, we have adapted the gwst code to compute the 
matrix elements of the screened Coulomb interaction in 
Eq. [31 In contrast with earlier work for semiconductors, 
the Auger matrix elements have been determined entirely 
from first principles. The calculations were performed at 
the experimental lattice constants reported in Ref. ITol. 
For additional technical details and convergence parame- 
ters we refer to Refs.|7| and^. We also note that a plane- 
wave cutoff of 17 Ha is required for the matrix elements 
such that the Auger rate is converged to within 10%. 

For both eeh and hhe Auger processes we distin- 
guish between intra- and mter-band events, depending 
on whether the final state of the scattered carrier lies in 
the same band as the initial state (see insets of Fig. [T]). 
For typical concentrations of injected carriers (10^^-10^° 
cm~^) holes and electrons are confined to a small region 
of the Brillouin zone around the F point. Assuming equal 
electron and hole concentrations the maximum momen- 
tum transfer to the Auger electron or hole that is scat- 
tered into the final state is 'ikp (at T = OK), where kp 
is the Fermi momentum). Since energy and momentum 
have to be conserved in the Auger process, intraband 
Auger events are therefore only likely to occur for alloy 
compositions with small band gaps, as we demonstrate 
later. For larger band gaps intraband processes are neg- 
ligible, and the possibility of interband Auger processes 
then depends on whether other valence or conduction 
bands are available into which holes or electrons can be 
scattered while energy is conserved. Our accurate band 
structures allow us to determine that possibility. 



Focusing first on eeh processes, we note that the low- 
est conduction band in all three binary wurtzite phases is 
nondegenerate and almost spherically symmetric around 
the F-point. The next higher conduction band is situated 
at an energy A above the CBM at F ranging from 2.5 to 
3.1 eV (see Fig. [2] and Table [I]). This strongly suggests 
that interband eeh Auger processes that excite an elec- 
tron into this second conduction may occur for InGaN 
alloys with matching band gaps. This second conduction 
band was not included in a recent k.p study by Hader et 
ai, leading them to conclude that direct Auger losses are 
negligible in InGaN quantum wells piS, 

The formidable computational challenges involved in 
calculating Auger rates were discussed by Laks et ali^ 
The calculation of the Auger recombination rate itself 
[Eq. IH)] scales with the fourth power of the number of 
k points, rendering explicit integration over a grid of 
k points prohibitively expensive. We therefore tackled 
this multidimensional integration with a Monte Carlo ap- 
proach. We compute statistical averages over 40,000,000 
Monte Carlo steps, chosen so that the error of the mean 
is always at least one order of magnitude lower than the 
value of the mean. Because of the expense involved in 
a first-principles calculation of Auger matrix elements 
[Eq. (121)], we confined the calculations to a finite num- 
ber of points in the twelve-dimensional fc-space (10,000 
elements in total), and employed a linear-interpolation 
scheme for k points off the mesh. This approach is jus- 
tified because we found that the matrix elements vary 
only weakly over the small volume of the Brillouin zone 
involved in Auger recombination. Accurate interpola- 
tions of the OEPx(cLDA)-|-GoWo band structures were 
obtained using an anisotropic effective-mass modelj^ 

Our present calculations are aimed at examining Auger 
rates for a wide range of InGaN band gaps. Explicit eval- 
uation of band structures and wave functions for each 
of the corresponding alloy compositions would be pro- 
hibitive. We performed calculations entirely from first- 
principles for eeh processes in pure GaN, and modeled 
different alloy compositions by applying a "scissor shift" 
to the band gap (i.e., rigidly shifting all conduction bands 
relative to the valence bands). While this choice is not 
optimal for calculations of alloys with large In concentra- 
tions, we note that In concentrations in current optoelec- 
tronic devices do not exceed 15% and thus an extrapola- 
tion of band parameters and matrix elements from GaN 
is the most sensible choice. With each Eg value we asso- 
ciate an alloy composition by using a bowing parameter 
b — 2.5eVji2, and then use linear interpolation to obtain 
A based on the values listed in Table HI This approach is 
justified by the similarity of the wave functions and band 
structure for corresponding bands of GaN and InN. 

The Auger coefficient, defined as C — R/n^, is re- 
ported in Fig. [1] as a function of the fundamental band 
gap for T — 300 K and a carrier density of 1 x 10^^ cm~'^. 
For the purposes of the figure, we populated only the 
light-hole band with holes. Holes in other valence bands 
lead to a similar magnitude of the Auger coefhcient. In 
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6.47 


3.24 


0.69 


6.53 3.07 


0.53 
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2.65 


2.50 


3.12 


8.62 9.07 


10.12 


Ae 


-0.295 


0.034 


0.066 







TABLE I: Band gap [Eg), gap between the first and the sec- 
ond conduction band at V (A) and crystal-field splitting (Ac) 
calculated with the OEPx(cLDA)-|-GoWo approach for AIN, 
GaN and InN. 



our calculations, which use Fermi-Dirac statistics, the 
Auger coefficient is largely independent of n for densi- 
ties up to 1 X lO^^cm"'^. We note that the interband 
process has a weak dependence on the choice of bowing 
parameter due to the method of choosing A. For h rang- 
ing from 1 to 4 eV, the center of the recombination peak 
varies by less than 0.1 eV. 

Three distinctly different regimes emerge. For large 
band gaps, intraband Auger recombination is negligible 
and interband recombination dominates. For small band 
gaps (< 1.0 eV) the situation is reversed. In an inter- 
mediate energy region, both eeh recombination processes 
are irrelevant. The two Auger regimes exhibit a differ- 
ent characteristic behavior as a function of band gap: 
while interband recombination peaks when A and Eg are 
in resonance, intraband recombination monotonically in- 
creases with decreasing band gap. hhe processes are not 
included here, but inspection of the band structures in 
Fig. [5] allows us to conclude that they do not contribute 
at Eg values above 2.5 eV. 

The Auger coefficients computed with our model fall 
within the (very wide) experimentally reported range of 
1x10-3"* to 5x10-2* cm^s-i (Ref.S). In particular, they 
agree well with the values reported by Shen et al^ based 
on their optical pumping experiments. Their reported 
Auger coefficient ranged from (1.4 - 2.0) xlO-^^ cm^s-* 



and increased with In concentrations rising from 9 to 
15%. Our calculated values in Fig. [1] are in qualita- 
tive and quantitative agreement with these observations, 
strengthening the case that the observed losses are indeed 
due to Auger recombination. Interestingly, Fig. [T] shows 
that the Auger coefficient continues to increase when the 
band gap is lowered down to 2.5 eV; i.e., for longer wave- 
lengths. This indicates that rising Auger losses could well 
play an important role in the "green gap" problem, the 
well known decrease in efficiency of InGaN light emitters 
at longer emission wavelengths. 

Based on our results we can examine possible strategies 
for reducing Auger losses in the blue/green region of the 
spectrum. (1) Using the zinc-blende phase, in which the 
second conduction band occurs at much higher energies 
(see Table IJ); however, growth of high quality phase-pure 
zinc-blende nitrides has proven very difficult. (2) Strain 
engineering the band structure to move the second con- 
duction band away from the resonance; however, both 
the first and second conduction band are nondegenerate 
and thus sensitive only to hydrostatic strain, and our cal- 
culated deformation potential for the energy difference A 
is quite small. (3) Tuning InGaAlN alloy compositions. 
Again, this is likely to be fruitless because the value of A 
is similar for all three nitrides (Table IJ). 

In summary, we have presented first-principles eval- 
uations of Auger recombination rates in InGaN alloys, 
showing the presence of a resonance that leads to in- 
creasing values of the Auger coefficient for wavelengths 
ranging from blue to green. The calculated values are in 
good agreement with experiment and confirm that Auger 
recombination is a key loss mechanism in nitride light 
emitters. 
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